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Abstract 

This paper considers a disturbance attenuation problem for a linear discrete time 
invariant system under random disturbances with imprecisely known distributions. The 
statistical uncertainty is measured in terms of the mean anisotropy functional. The dis- 
turbance attenuation capabilities of the system are quantified by the anisotropic norm 
which is a stochastic counterpart of the Tioo norm. The designed anisotropic subopti- 
mal controller generally is a dynamic fixed-order output-feedback compensator which is 
required to stabilize the closed-loop system and keep its anisotropic norm below a pre- 
scribed threshold value. Rather than resulting in a unique controller, the suboptimal 
design procedure yields a family of controllers, thus providing freedom to impose some 
additional performance specifications on the closed-loop system. The general fixed-order 
synthesis procedure employs solving a convex inequality on the determinant of a posi- 
tive definite matrix and two linear matrix inequalities in reciprocal matrices which make 
the general optimization problem nonconvex. By applying the known standard convexi- 
fication procedures it is shown that the resulting optimization problem is convex for the 
full- information state- feedback, output-feedback full-order controllers, and static output- 
feedback controller for some specific classes of plants defined by certain structural proper- 
ties. In the convex cases, the anisotropic 7-optimal controllers are obtained by minimizing 
the squared norm threshold value subject to convex constraints. In a sense, the anisotropic 
controller seems to offer a promising and flexible trade-off between %2 and Hoc controllers 
which are its limiting cases. In comparison with the state-space solution to the anisotropic 
optimal controller synthesis problem presented before which results in a unique full-order 
estimator-based controller defined by a complex system of cross-coupled nonlinear ma- 
trix algebraic equations, the proposed optimization-based approach is novel and does not 
require developing specific homotopy-like computational algorithms. 
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1 Introduction 



The stochastic uncertainty of random disturbances regarded as a discrepancy between an in- 
exactly known probability distribution of a real-world noise and its nominal model can signifi- 
cantly degrade the designed performance of a control system if the applied controller synthesis 
procedure relies upon a specific probability law of the disturbance and the assumption that it is 
known precisely. Such situations can also result from the inherent variability of the conditions 
of the control system operational environment. So, the I-L2 and Hoo controllers are efficient in 
full only if the basic hypotheses on the nature of external disturbances are met closely enough. 
It is known that the H2 (or LQG) controller may perform poorly if the input disturbance is 
a strongly correlated noise [TJ, while the controller designed for the deterministic worst 
case [2] demonstrates excessive conservatism if the external disturbance is white or weakly 
correlated random signal. 

One of the first ideas aimed at overcoming the lack of performance of the LQG controller 
in the case when the external disturbance is not the Gaussian white noise arose in work [3] 
devoted to some modification of the performance criterion. This idea gave rise to development 
of the whole class of problems in the control theory called the risk sensitivity problems [U [5] . 

The ideas of deriving controller which combines the positive features of LQG (H2) and Hoo 
controllers (i.e. minimizes the quadratic cost sufficiently good and is robust enough) appeared 
in the beginning of 1990's. In particular, one can distinguish an approach concerned with 
minimization of H.2 norm of the closed-loop system under constraints on its norm [6J and 
approach related to minimization of "Hoc entropy functional under constraints on the closed-loop 
Hoo norm [7j. 

As is shown in [8j , the problem of synthesis of a controller which minimizes the entropy 
functional is equivalent in a sense to the problem of optimal risk-sensitive (LEQG) controller 
synthesis. A lot of papers are devoted to the problems concerned with minimization of the H,^ 
entropy functional (see e.g. [9j-[13j). 

The ideas of the mixed T^/'Hoo control first introduced in [6j were extended in [THEE] based 
on splitting the external disturbance into signals with bounded spectrum and bounded power 
and using the multi-objective %2/'Hoo performance criterion. A solution to the stochastic mixed 
^2/^00 control problem for the discrete-time systems is given in |16j . 

All of the works mentioned above exploit the techniques based on solving certain (sometimes 
cross-coupled) Riccati equations. In [T7j the mixed H2/H00 problem was considered in terms of 
algebraic Riccati inequalities rather than equations and solved by means of convex optimization. 
Since then the efficient interior-point algorithms for solving convex optimization problems had 
been developed [IH] [20] , convex optimization has become a standard strategy for control system 
analysis and synthesis. The linear matrix inequalities have proved to be a powerful formulation 
and design technique for a variety of linear problems [21j . After the Hoo controller synthesis 
problem had been solved via LMI [23l [22] . the semidefinite programming was successfully 
applied to developing effective solutions to multi-objective ^2/^00 control problems [2"3]-[34|. 
A detailed survey of these extensive results is far beyond the topic of this paper and may be 
presented elsewhere. 

An approach to attenuation of uncertain stochastic disturbances based on minimax control 
was proposed in the middle of 1990 's [35] [37] and extended later to the MIMO systems and 
synthesis of structured controllers via LMI in [35]. Instead of exact knowledge of the distur- 
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bailee's covariance coefficients, it is only required that the covariance coefficients belong to an a 
priori known set. The designed controller minimizes the worst possible asymptotic output vari- 
ance for all these disturbances. The considered problem is intermediate between the extreme 
%2 and %oo design scenarios and reduces to a robust control problem with uncertainty in the 
external disturbance signal [38J. 

At the same time, another promising stochastic minimax alternative had emerged from ideas 
of I.G. Vladimirov who originally developed the anisotropy-based theory of robust stochastic 
control presented in a series of papers [12] . In the view of this approach, the robustness 
in stochastic control is achieved by explicitly incorporating different scenarios of the noise 
distribution into a single performance index to be optimized; the statistical uncertainty is 
measured in entropy theoretic terms, and the robust performance index can be chosen so as to 
quantify the worst-case disturbance attenuation capabilities of the system. The main concepts 
of the anisotropy-based approach to robust stochastic control are the anisotropy of a random 
vector and anisotropic norm of a system. 

The anisotropy functional introduced by I.G. Vladimirov is an entropy theoretic measure 
of the deviation of a probability distribution in Euclidean space from Gaussian distributions 
with zero mean and scalar covariance matrices. The mean anisotropy of a stationary random 
sequence is defined as the anisotropy production rate per time step for long segments of the 
sequence. In application to random disturbances, the mean anisotropy describes the amount of 
statistical uncertainty which is understood as the discrepancy between the imprecisely known 
actual noise distribution and the family of nominal models which consider the disturbance to 
be a stationary Gaussian white noise sequence with a scalar covariance matrix [12| H3"] . 

Another fundamental concept of I.G. Vladimir ov's theory is the a-anisotropic norm of a 
linear discrete time invariant (LDTI) system which quantifies the disturbance attenuation ca- 
pabilities by the largest ratio of the power norm of the system output to that of the input 
provided that the mean anisotropy of the input disturbance does not exceed a given nonneg- 
ative level a (HI H3]. A generalization of the anisotropy-based robust performance analysis to 
finite horizon time varying systems is developed in [H]. 

In the context of robust stochastic control design aimed at suppressing the potentially 
harmful effects of statistical uncertainty, the anisotropy-based approach offers an important 
alternative to those control design procedures that rely on a precisely known specific probabil- 
ity law of the disturbance and the assumption that it is known precisely. Minimization of the 
anisotropic norm of the closed-loop system as a performance criterion leads to internally stabi- 
lizing dynamic output-feedback controllers that are less conservative than the Hoc, controllers 
and more efficient for attenuating the correlated disturbances than the I-L2 controllers [43J. 
A state-space solution to the anisotropic optimal control problem derived by I.G. Vladimirov 
in jl5] involves the solution of three cross-coupled algebraic Riccati equations, an algebraic Lya- 
punov equation and an equation on the determinant of a related matrix. The resulted optimal 
full-order estimator-based (central) controller is a unique one. An extension of these results 
to the systems with parametric uncertainties was given in E] • But solving these complex 
systems of equations requires special developing of homotopy-like numerical algorithms [35] • Be- 
sides, the applied equation-based synthesis procedure is not aimed at the synthesis of reduced- 
or fixed-order (decentralized, structured, multi-objective) controllers which still remains open. 
Moreover, although the ideas of entropy-constrained induced norms and associated stochastic 
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minimax find further development in the control literature [49] . the anisotropy-based theory 
of stochastic robust control remains largely unnoticed. One of the reasons seems to be hard 
numerical tractability of the analysis and synthesis problems as well as a lack of additional 
degrees of freedom in the controller synthesis procedure. 

The anisotropic suboptimal controller design is a natural extension of the approach proposed 
by I.G. Vladimirov in [15]. Instead of minimizing the anisotropic norm of the closed-loop system, 
a suboptimal controller is only required to keep it below a given threshold value. Rather than 
resulting in a unique controller, the suboptimal synthesis yields a family of controllers, thus pro- 
viding freedom to impose some additional specifications on the closed-loop system. One of such 
specifications, for example, may be a particular pole placement to achieve desirable transient 
performance. Getting a solution to the anisotropic suboptimal controller synthesis problem 
requires a state-space criterion to verify whether the anisotropic norm of a system does not 
exceed a given value. An Anisotropic Norm Bounded Real Lemma (ANBRL) as a stochastic 
counterpart of the well-known T-L^ norm Bounded Real Lemma for LDTI systems under statis- 
tically uncertain stationary Gaussian random disturbances with limited mean anisotropy was 
presented in [50]. The resulting criterion has the form of an inequality on the determinant of a 
matrix associated with an algebraic Riccati equation which depends on a scalar parameter. A 
similar criterion for linear discrete time varying systems involving a time-dependent inequality 
and difference Riccati equation is derived in [5Tj . Recently, a sufficient strict version of ANBRL 
was introduced in [52j [53] in form of a convex feasibility problem employing a strict inequality 
in the determinant of a positive-definite matrix and a related LMI. Moreover, the determinant 
constraint turns out to depend linearly on the squared threshold value, thus allowing to min- 
imize it directly subject to the convex constraints and compute the a-anisotropic norm of a 
LDTI system as a solution to the convex optimization problem [53]. The developed analysis 
procedure is numerically attractive and easily realizable by means of available convex optimiza- 
tion software [34], [55] . This paper is aimed at application of the powerful technique of convex 
optimization to synthesis of the anisotropic suboptimal and 7-optimal controllers generally of 
fixed order. The anisotropic controller seems to offer a promising and flexible trade-off between 
%2 and "Hoc controllers. In comparison with the state-space solution to anisotropic optimal 
controller synthesis problem derived before in [33], the proposed optimization-based approach 
is novel and does not require developing specific homotopy-like computational algorithms |48j. 

The structure of the paper is as follows. In Section [2] we give the statement of the general 
problem of synthesis of the fixed-order anisotropic suboptimal controller. In Section [3] we in- 
troduce a solution to the general fixed-order synthesis problem and consider three important 
design cases: static state-feedback gain for full-information case, dynamic output-feedback con- 
troller, and static output-feedback gain. Section [31 provides a number of illustrative numerical 
examples. Concluding remarks are given in Section [3J 

1.1 Notation 

The set of reals is denoted by R, the set of real (n x m)- matrices is denoted by R nxm . For a 
complex matrix M = [m^], M* denotes the Hermitian conjugate of the matrix: M* := [m*J. 
For a real matrix M = [my], M T denotes the transpose of the matrix: M T := [m^]. For real 
symmetric matrices, M y N stands for positive definiteness of M — N. In block symmetric 
matrices, symbol * replaces blocks that are readily inferred by symmetry. The spectral radius of 
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a matrix M is denoted by p(M) := max/j |Afc(M)|, where Afc(M) is k-th eigenvalue of the matrix 
M. The maximum singular value of a complex matrix M is denoted by a(M) := ^/A max (M*M). 
7 n denotes (n x ra) identity matrix, nxm denotes zero (n x m) matrix. The dimensions of zero 
matrices, where they can be understood from the context, will be omitted for the sake of brevity. 

The angular boundary value of a transfer function F(z) analytic in the unit disc of the 
complex plane \z\ < 1 is denoted by 

F(u) := lim F(re itJ ). 

r— 

T-L^ m denotes the Hardy space of (p x m)-matrix-valued transfer functions F(z) of a complex 
variable z which are analytic in the unit disc \z\ < 1 and have bounded %2 norm 

1/2 

/if-- \ 

^2 



— / tr(F(u)F*(u))du 



W* m denotes the Hardy space of (p x m)-matrix-valued transfer functions F(z) of a complex 
variable z which are analytic in the unit disc \z\ < 1 and have bounded "Hoc norm 



sup a(F(z)) = ess sup g(F(uj)). 



2 Problem statement 

Consider a LDTI plant P(z) with n^-dimensional internal state X, m w -dimensional disturbance 
input W, m u -dimensional control input U, p 2 -dimensional controlled output Z, and ^-dimen- 
sional measured output Y. All these signals are double-sided discrete-time sequences related to 
each other by the equations 
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(1) 



where all matrices are assumed to be of appropriate dimensions and p z ^ m w ; (A,B U ) and 
(A, C y ) are assumed to be stabilizable and detectable. 

The only prior information on the probability distribution of the disturbance sequence 
W = (wk)-oo<k<+oo is as follows. It is assumed that W is a stationary sequence of random 
vectors Wk with zero mean Ewk = 0, unknown covariance matrix EwkiuJ = Tiw >- 0, and 
Gaussian PDF 

1 



p(w k ) := (2tt) 



-mi/2 



(detS 



w) 



-V2, 



cxp --\\w k \ 



where 



w^Ti^Wk and E denotes the expectation. It is also assumed that the mean 
anisotropy of the sequence W is bounded by a nonnegative parameter a. The latter means that 
W can be produced from m w -dimensional Gaussian white noise V = (v k)-oo<k<+<x> with zero 
mean Ei> fe = and scalar covariance matrix Ew fc f J = XI mi , A > 0, by an unknown stable LTI 
shaping filter G(z) in the family 

G a := {GeH? xm : A(G)^a}, 
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where 



A(G) = ~ I lndet U^-G{u)G*{u) ) du; 



is the mean anisotropy functional [I21 

We are generally interested in finding a fixed-order dynamic output-feedback controller in 
general compensator form 



K(z): 
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-oo < k < +00, 



(2) 



with rig-dimensional internal state E = (£fc)_oo<fc<+oo to ensure stability of the closed- loop 
system (Figure Q]) and guarantee some designed level of the external disturbance attenuation 
performance. 



Z< 

Y 



W 



u 



K 



Figure 1: Closed- loop system 

Let Tzw( z ) denote the closed-loop transfer function from W to Z. Recall that the a- 
anisotropic norm of a transfer function Tzw(z) £ H^* 171 ™ quantifies the disturbance attenuation 
capabilities of the respective closed-loop system by the largest ratio of the power norm of the 
system output to that of the input provided that the mean anisotropy of the input disturbance 
does not exceed the level a [321 S3] : 
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(3) 



Moreover, it is known from 



that the a-anisotropic norm of a given system F G "H^ 



is a nondecreasing continuous function of the mean anisotropy level a which satisfies 

1 
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zw\\a ^ hm |T^|| 

a— >+oo 



T; 
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(4) 



These relations show that the %2 and "Hoo norms are the limiting cases of the a-anisotropic 
norm as a — > 0, +00, respectively. 

The statement of the general problem of synthesis of the fixed-order anisotropic suboptimal 
controller is as follows. 

Problem 1. Given a LDTI plant P described by (TJP, a mean anisotropy level a ^ of the 
external disturbance W , and some designed threshold value 7 > ; find a fixed-order LDTI 
output-feedback controller K defined by |||) which internally stabilizes the closed-loop system 
and ensures its a-anisotropic norm does not exceed the threshold 7, i.e. 



IT 
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(5) 
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3 Problem solution 



Here we introduce a solution to the general fixed-order synthesis problem and consider three 
important design cases, namely static state-feedback gain for full-information case, dynamic 
output-feedback controller, and static output-feedback gain. To solve the synthesis problem, 
we apply a state-space criterion to verify if the anisotropic norm of a system does not exceed a 
given threshold value. This criterion called the Strict Anisotropic Norm Bounded Real Lemma 
(SANBRL) was recently presented in [52], 153]. But to apply SANBRL to the synthesis problem 
we should recast it in slightly different form. 



3.1 Anisotropic norm bounded real lemma 

With the plant P and controller K defined as above, the closed-loop system admits the real- 
ization 



Tzw(z) : 
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-oo < k < +oo, 



(6) 



where Xk € K n , n = n x + n^. 

It is shown in J52J [53] that given a ^ 0, 7 > 0, the inequality holds true if there exists 
77 > 7 2 such that the inequality 



£ T $2 - £> T 2)))) 1/mu ' < 7 2 



holds for a real (n x n)-matrix $ = <3> T y satisfying LMI 

a t $a - $ + e T e yi T $s + c t d 



(7) 



-< 0. 



Note that the constraints described by the inequalities ([7]) and ([H]) are convex with respect to 
both variables r] and $. Indeed, the function — (det(-)) 1 / m '" of a positive definite (m w x m w )- 
matrix on the left-hand side of (JTj) is convex; see [T9l 156] . 

Being convex in both variables r] and $, the conditions (JTj), ([8]) of SANBRL j52j [53] are not 
directly applicable to solving the intended synthesis problem because of the cross-products of the 
unknown Lyapunov matrix $ and the closed-loop realization matrices (A, 25, C, D) depending 
affinely on the controller parameters, which also appear in ([7j). Moreover, just the inequality ([7]) 
does not allow for the well-known Projection Lemma [221 E3] to be applied to get rid of the 
controller realization matrices in the synthesis inequalities. 

To overcome this obstacle, let us first move the positive definite matrix T]I mw — S T $!B — D T D 
away from the determinant in ([7J by introducing a slack variable, real (m w x m^)-matrix 
^ = y such that 



7] - (det (e- 2 ^™-^)) 1 /^ < 7 a j ^ _< ^ 



(9) 



which is equivalent to (JTJ). Then, let us decouple the cross-products of $, H>, and D in Qj. For 
this purpose, the latter inequality in ([9]) can be rewritten as 
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where 



.$-1 




-< 0, which is equivalent to 
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by virtue of the Schur Theorem; see e.g. [21 
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where — $ ^ evidently. Then by the Schur Theorem the last inequality is equivalent to 
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To decouple the cross-products of C and D, let us represent the inequality ffTTj) as 
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where — I Vz -< clearly. Second application of the Schur Theorem to the above inequality gives 
the following formulation of SANBRL in reciprocal matrices. 

Lemma 1. LetTzw £ 'H^ 1 xrriw be a system with the state-space realization (T2|), where p(A) < 1. 
T/ien zts a- anisotropic norm (T3j] is strictly bounded by a given threshold 7 > 0, i.e. HT^f^ < 7 
if there exists r\ > 7 2 stzc/i £/ia£ £/ie inequality 



77- (det ( e - 2a /m^))l/m» <7 2 



:i2i 



/ioWs true /or some rea/ (m w x m w )-matrix \I/ = \I/ T >- and (n x n)-matrix $ = $ T >- 
satisfying inequalities 
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(13) 



(14) 



Thus, with the notation II := $ , verifying if the condition |||TzvHIL < 7 holds true reduces 
to finding a positive scalar 77 and two matrices $ >- 0, II >- 0, $11 = I n , satisfying the LMIs fll3p . 
(I14p under the convex constraint (j!2p or making sure of insolvability of this problem. For solving 
this nonconvex problem numerically, one can make use of known algorithms developed in |57j- 
[65] suitable for finding reciprocal matrices under convex constraints. 
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3.2 State-feedback controller 



Before to proceed to general synthesis Problem dj let us consider the full- information case, when 
the state vector can be measured precisely and the plant is described by the equations 
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(15) 



where as above all matrices are assumed to be of appropriate dimensions and p z ^ m w ; (A, B u ) 
is assumed to be stabilizable. 

Problem 2. Given a L DTI plant P described by / f75]) . a mean anisotropy level a ^ of the 
external disturbance W , and some designed threshold value 7 > ; find a static state-feedback 
controller 

u k = Kx k (16) 
which internally stabilizes the closed-loop system Tzw{ z ) with the state-space realization 



(17) 



and ensures its a-anisotropic norm does not exceed the threshold 7, i.e. the inequality holds. 

The following theorem gives sufficient conditions for the static state-feedback anisotropic 
suboptimal controller to exist. 

Theorem 1. Given a ^ 0, 7 > ; the state-feedback controller (T76j] stabilizing the closed-loop 
system ( fi7| ) (p(A + B U K) < 1 ) and ensuring (EJ) exists if the convex problem 

7] - (det [ e -^/m w ^y/m w < 
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(18) 
(19) 

(20) 

77 > 7 2 , * y 0, n >- (21) 

zs feasible with respect to the scalar variable 7], real (m w x m w )-matrix real (n x x n x )-matrix 
IT, end rea/ (m u x n x ) -matrix A. // the problem (f7%j ) -(ff7| ) feasible and the unknown variables 
have been found, then the state-feedback controller gain matrix is determined by K = AIT -1 . 
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Proof. Let a solution to the problem (fT8"l) - fl2~Tl) exist. Define $ 
An" 1 , the LMIs (fT3]). (|2Djl can be rewritten as 
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Pre- and post-multiplying the last inequality by blockdiag($, I mw , I nx , L Pz ) y yields 
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(23) 



Then, by Lemma HJ from (Tl8j) . (j22|) . (|23|) . (12T]) it follows that the controller gain matrix K is 
the solution to Problem [2] for the closed- loop realization (TTTj) , which completes the proof. □ 

Remark 1. Although it is not hard to prove that the synthesis inequalities lffi^) - i{2l\) and the 
conditions [W\) -§JJ$ of Lemma U\ are equivalent, we can only establish and prove sufficient 
existence conditions for the controller [To]) since the conditions of Lemma[J\ are only sufficient. 
This also concerns two further synthesis theorems. 

Corollary 1. The inequalities [W\) - ^21\) are not only convex in \1/ and affine with respect to 
IT and A, but also linear in r y 2 . Obviously, minimizing r y 2 under the convex constraints (TT$J)- 
( TJ5P ; we minimize 7 under the same constraints. With the notation 7 := 7 2 ; the conditions of 
Theorem^ allow to compute the minimal 7 via solving the convex optimization problem 



minimize 7 
over \l>, LT, A, r?, 7 satisfying ffity -[21 



(24) 



// the convex problem §2J$ is solvable, the state-feedback controller gain matrix is constructed 
just as in Theorem^ 

All anisotropic controllers obtained from solutions to optimization problems like (T24"]) will 
be referred to as anisotropic 7- optimal controllers. 

3.3 Fixed-order output-feedback controller design: 
convex constraints on reciprocal matrices 

Direct application of the sufficient conditions ffl~2T) — (TT4|) of Lemma [1] to the closed- loop realiza- 
tion 
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yields the following corollary on the straightforward solution to general Problem [TJ 

Corollary 2. Given a ^ ; 7 > ; a dynamic output-feedback controller K of order defined 
by (0) solving Problem^ exists if the inequalities 
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are feasible with respect to the scalar variable r\, real (m w Xm w ) -matrix ^ , matrices A c 6 
£? c G R n « xp «, C c E ]R m « xn e ; £) c g ]R m " xp f and two reciprocal (n x n) -matrices $, II suc/i that 



where n = n x + is the closed-loop system order. 



(30) 



Thus, the problem of finding the realization matrices (A c , B c , C c , D c ) of the fixed-order 
output-feedback dynamic controller (j2J) solving Problem [1] leads to solving the problem ([26]) - 
fl30|) or making sure of its insolvability. The problem f[2"6"]) - f[30"]) is nonconvex because of the 
condition fl30l) . Although application of the known algorithms of [57] [65] can leads to a suc- 
cessful solution of the problem (j2"6"l) - (j501) . it should be kept in mind that any of them can 
converge to local minima. Nevertheless, the full-order controller synthesis allows for a quite 
standard convexification procedure which is considered below to be applied. 



3.4 Full-order output-feedback controller 

For full-order design (n x = rig) one can effectively apply the well-known linearizing change of 
variables presented in [66J and used in [26] in the multi-objective control framework. From the 
block partitioning in (I2"9"j) and the condition (|3"Uj) it follows that 



which leads to 



with the notation 



$i := 



$rii = $i, 

In. $11 

o $? 2 



" n n " 




" In 










n$i = iii 
nn 



ni ■-- 



uj 2 



It can be easily shown by direct calculation that 



(31) 



(32) 



$^n$i = iTf 



In x $11 



(33) 



11 



The key linearizing change of the controller variables is defined as follows [66] 

a c := $i 2 A c nT +$ 12J B c c 2/ n 11 + $ llJ B u c c n? 2 + $ 11 (A + J B^ c c 2/ )n 11 , 



2), 



+ D c c y u 

A,. 



(34) 
(35) 
(36) 
(37) 



The new variables A c , S c , C c , 2) c have dimensions n x x n x , n x x p y , m u x n x , and m u x 
respectively, even if n x ^ n^. It is noted in [26] that if $12 and IT 12 have full row rank and if 
A c , S c , C c , T> c , riii, and $n are known, one can always find the controller matrices A c , I? c , C c , 
D c satisfying (I3"4"I) - (13"T|) . If the matrices $ i2 and Il 12 are square {n x = n^) and invertible, then 
A c , B c , C c , and D c are unique, i.e. for full-order design, when one can always assume that $i 2 
and Ili 2 have full row rank, the mapping defined by fl34l) - fl37j) is bijective. More details can be 
found in [66], [26]. 

The solution to Problem [T] in the full-order design case is given by 

Theorem 2. Given a ^ ; 7 > ; a dynamic output-feedback controller K of full order = n x 
defined by (TJ|) solving Problem\J\ exists if the convex problem 



T) - (det ( e - 2a /«^))iAn„ < 7 2 , 



-nn 



AU n + B U G C 
A c 

C z Tlu 4- D ZU G C 



B w -)- B u T) c Dy W 
®i\B w + ( B c D yw 
B zw B ZU D c Dy W 



* 


A + B u V c C y 

+ T> C Cy 

C z + D zu T) c Cy 



* 
* 

-$11 





B„ 



* 

-nn 



B U D c By W 



®uB w + 'BcBy-w 
B zw 



B zi/D C Dyyj 



* 
* 

Pz 

* 
* 
* 

-nn 





77 > 7 2 , n n >- 0, $11 y 0, 



nn 4, 
4, $11 



-<o, 



* 



>- 



(3* 



(39) 



* 
* 
* 



-<0, 



(40) 



(41) 



is feasible with respect to the scalar variable t], real (m w x m w )-matrix *Sf, matrices A c G 



pn x xn x 



\n x xp y 



i>m u xri; 



x , £> c G ]R m " x Pf and two rea/ (n x x n x )-matrices n n , $ n . // £/ie 



problem (f^)-(^ip zs feasible and the unknown variables have been found, then the controller 
matrices A c G R n * xn * ; S c e R n * x J>« ; C c e l 1 ^^^ £> c G ]R m " x ^ are uniquely defined by 



B c 
A- 



2) r 



(e, 



D c c , 2/ n 11 )n 12 T , 



$uB u D c 



$r 2 (-A c - $i 2 5 c c y n n - ^n^Cen^ - $11 (A + s^c^niOnr/ 

and determined from finding two nonsingular (n x x n x )-matrices IT 12; $ i2 t/iat satisfy 

ni 2 ^ 2 = / n!c _nii$ii. 



(42) 
(43) 
(44) 
(45) 

(46) 
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Proof. Let a solution to (I5B]) - MI]) exist. From Ij32|) - (IB7]1 and (J25]) it follows that 



Ann + 5 u e c A + B u T) c C y 
A c $ n 7L + ScCj, 



$-f>uii, 



&iiB w + ¥> c D yw 



[ C z U.ii + D ZU G C C z + D zu T> c C y ] — CIL, 



nn 4, 
4, $n 



where $ and II are defined by (129]) and satisfy (130]) with = n x . Substitution of the above 
identities to the inequalities ( 139|) . (H0|) yields 



D 



-$^n$a 







^0, 



-n^$ni 







n^e T 




















cil 


D 





Vz 



^0. (47) 



Performing a congruence transformation with 

blockdiag(J mi „, $7 T , 7 P J, blockdiag(Il7 



; ^J7Hu J ^1 5 -^Pz 



on the inequalities (1171) . respectively, leads to 
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-n 














. 




Then, by Lemma (H from 


(EHD, 


(BHD, d 


29}, (J30D 



-$ 







e T 







B T 








-n 





e 










-< 0. 



(48) 



internally stable and its a-anisotropic norm does not exceed the designed threshold 7, i.e. the in- 
equality ([5]) holds. The procedure of reconstruction of the controller realization (A c , B c , C c , D c ) 
from the solution variables (A c , S c , C c , D c ) by gB}, (pj) - (j!5)l is quite standard jSS], |2S]. □ 



Corollary 3. As the inequalities (E2J)-([PP are also linear in 7 := 7 2 ; the conditions of Theo- 
rem^ allow to compute the minimal^ via solving the convex optimization problem 



minimize 7 

over $ u , nn, A;,B C , C c , 2^,77, 7 satisfying (Egj-CT. 



(49) 



7/ i/ie convex problem ( f7P| ) solvable, the controller matrices are constructed just as in Theo- 
rem^ 

It is stressed in [26] that the applied synthesis procedure does not introduce any conser- 
vatism, if the analysis result does not involve any. 

The results of Theorem [2] make possible application of the anisotropic norm as a closed- loop 
performance specification or objective for specific closed-loop channels in the multi-objective 
control problems based on a common Lyapunov functions (2H] together with other performance 
specifications and objectives that can be captured in the LMI framework. 
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3.5 Static output-feedback controller 

Let us now consider the special and very important case of static output-feedback controller 

u k = Ky k . (50) 

Problem 3. Given LDTI plant P described by (TJP, a mean anisotropy level a ^ of the 
external disturbance W , and some designed threshold value 7 > ; find the static output-feedback 
controller |37|) which internally stabilizes the closed-loop system Tzw{ z ) with the state-space 
realization 

(51) 
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3 " 




A + B u KC y 


B w -\- B U K Dy W 
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D 




C z + D zu KC y 


D zw ~\~ D ZU K Dy W 



and ensures its a-anisotropic norm does not exceed the threshold 7, i.e. 

\\Tzw\\ a < 7- 



(52) 



Direct application of the sufficient conditions ( ll2p - (fH"|) of Lemma [T] to the closed-loop 
realization (1511) yields the following corollary on the straightforward solution to Problem [3J 

Corollary 4. Given a ^ 0, 7 > 0, the static output-feedback controller / f50j) solving Problem^ 
exists if the inequalities 

(53) 
(54) 





71- (det ( e - 2a /™^))i/m u 
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2 

? 






* 






B w + B u KD yw — 11 
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Pz J 
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-n 


* 


C z -\- D ZU K Cy D zw -\- D ZU K Dy W 








-< 0, 



V > 7 



^^0, $^0, n^o 



(55) 



(56) 



are feasible with respect to the scalar variable r], real (m w x m w )-matrix real (m u xp y )-matrix 
K , and two reciprocal real (n x x n x )-matrices IT such that 



(57) 



So, the problem of finding the output-feedback gain matrix K solving Problem [3] leads 
to solving the problem (J5"3"]) - (|5"T]) or making sure of its insolvability. The inequalities (J5"3"]l - 
(|57p derived from the straightforward application of Lemma [T] are not convex because of the 
condition f|5T|) . One can try to solve this general problem by the algorithms of suitable 
for finding reciprocal matrices under convex constraints. 

However, the specific linearizing change of variables presented in [30] can make the result- 
ing optimization problem convex for a specific class of plants defined by a certain structural 
property. Namely, suppose that the transfer function of the plant (pQ) from the control input to 
measured output vanishes, i.e. 



T yu (z) := C y (zl - A)~ B. 



0. 



(5* 
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For the stabilizable and detectable plant (pQ), if (|58|) holds, then there exists a similarity trans- 
formation T such that 



' TAT -l 


TB W 






D zw 


D zu 


CyT^ 


Dy w 








A 12 


B W1 


B Ul 





A22 


B W2 





c zi 




D zw 


D zu 





Cy 2 


Dy w 






(59) 



where (A u ,B Ul ) is controllable, (A n ,C y2 ) is observable, and the matrix A 2 2 is stable [30J; see 
also [67] . The representation (|59i) implies that the closed-loop system realization after static 
output feedback becomes 



(60) 
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•B " 
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'A n 


A12 


+ B Ul KC y2 


B Wl 


+ B Ul KD yw 







A22 




B W2 
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+ D zu KCy 2 
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~\~ D ZU K D yw 



The Lyapunov matrix $ in the inequalities f fl3|) . ( Tb4"j) of Lemma [T] is partitioned according 
to the representation of A in fl6D|) as 



$ n $12 



>- 0. 



(61) 



The key linearizing change of variables is defined in 
P :-- 



as 
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s ' 




_ s T 
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12 



- $ 12$11 $22 - $i2$n $ 



12 



(62) 



" $11 


$12 " 




r q- 1 


-Q-'S 


.*£ 


$22 . 




_ S T Q- 1 


R-S T Q- 1 S 



It is noted in [30] that the transformation (1621) maps the set of all positive definite matrices into 
the set of all matrices with positive definite diagonal blocks; this map is bijective; its inverse is 
given by 

r A rT» 1 T /O-l /O-l C 1 

(63) 



(64) 
(65) 



The transformation fl62l) is motivated by the factorization [30] 
with 



Q 



P2 := 



/ -S 
i? 



Theorem 3. Suppose that the plant P described by (TJP is such that T yu (z) = 0, i.e. ( T55j) holds. 
Given a ^ 0, 7 > ; a static output-feedback controller defined by solving Problem^ exists 
if the convex problem 

(66) 



r? - (det ( e - 2 »K|r))iK < 7 2 ) 
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(67) 
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7] > 7 2 , * >- 0, 


QyO, RyQ 









-<0, 



(68) 
(69) 



is feasible with respect to the scalar variable r], real (m w x m w )-matrix , controller gain matrix 
K and real matrices Q, R, and S. 



Proof. Let a solution to flBBjl — exist. Then from (j55jl. (jBtljl . <IS5I) it follows that 

= Pi$Pi T , 



Q 
R 



A U Q A n S - SA 22 + A 12 + B U1 KC V2 
RA 22 
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RB W2 

[ C Zl Q C Zl S + C Z2 + D zu KCy 2 ] 



Pi$B, 



(70) 

(71) 

(72) 
(73) 



Substituting the identities (j7D]l - (j75]l to the LMIs HID, (JSHD, we have 



Pi$B -Pi^pT 
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-I 



Vz 



-<0, 



-p^pt o p^pT pe T 

- V I mw S T <&PT D T 
Pi^P^ Pi$S -Pi^pT 



GP 1 



T 
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-< 0. (74) 



Performing a congruence transformation with 

blockdiag(/ mu ,, Pf 1 , I p J, blockdiag(Pf \ 1^, Pf \ J 3 
where Pi is defined by (|65i) . on the inequalities (JHJ), respectively, yields 
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Pz J 










e T 
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-Ip. 



-< 0. 



(75) 



Pre- and post-multiplying the inequalities (f75|) by 

blockdiag(/ mm , J P J >- 0, blockdiag^, I mw , $~\ J P J >- 0, 
respectively, we have 



* - r?/ mw B T 2) T 
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(76) 
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Then, by Lemma (TJ from (1661) . (T76|) . (|69|) . (IBTj) it follows that the controller gain matrix K is 
the solution to Problem [3] for the plant (IBUj) and the closed- loop system ( loTij) . which completes 
the proof. □ 

Corollary 5. The convex constraints (G2|)-(ES)) are also linear in r y 2 . With the notation 7 := j 2 , 
the conditions of Theorem [3] allow for 7 to be minimized via solving the convex optimization 
problem 

minimize 7 
over Q, R, S, K, 77, 7 satisfying / Tff6)) -/ fffpp . 



(77) 



The controller gain matrix K enters the synthesis LMIs (167|) . (1681) directly. It is noted in [30] 
that this allows for some structural requirements on this controller gain to be incorporated 
making possible even the synthesis of decentralized controllers (with block-diagonal K) via 
convex optimization. 

The results of Theorem [3] make possible application of the anisotropic norm as a closed-loop 
performance specification or objective for specific closed-loop channels in the multi-objective 
control problems with LMI specifications considered in [30] . 

It should be also noted that in general case, when the structural property fl58|) does not hold, 
one can follow the way of [SU] and make use of the Youla-Kucera parametrization of stabilizing 
controller [5B1 ES] to parametrize aflinely the closed-loop system, enforce the said property, 
and bring the closed- loop realization to the form (159]) . Then the synthesis of the anisotropic 
controller can be treated as finding the Youla parameter that enters the closed-loop system 
aflinely by applying the results of Theorem [3] and Corollary [5j 

Besides the class of systems which satisfy the structural property (158 p , there are two partic- 
ular cases of the system's structure which allow for the static output-feedback design problem to 
lead to some convex optimization problem by applying a nonsingular state coordinate transfor- 
mation and introducing structured slack variables just as it was done for synthesis problem 
in [70]. These cases are the so called singular control and filtering problems. 

Let us first consider the singular control problem when the matrix D zu of the plant ([T]) is 
zero and the matrix B u is of full column rank. Then there exists a nonsingular state coordinate 
transformation matrix T u such that [70] 



(78) 



B u ■ T U B U 



Under this transformation, the plant realization matrices become 

A := T U AT U , B w := T U B W , C z := C Z T U , C y := C y T u . 



(79) 



Theorem 4. Suppose that the plant P described by ([!]) is such that D zu = and ranki? M = m u . 
Given a ^ ; 7 > ; a static output-feedback controller defined by [5D\) solving Problem^ for 
the closed-loop realization 



(80) 
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exists if the convex problem 



T) - (det (e^/™^)) 1 /™™ < 7 2 , 
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(82) 

(83) 
(84) 



?7 > 7% ^ ^ 0, $ >- 0, 

where A, B w , C z , C y are defined by (7S|) ; feasible with respect to the scalar variable 7], real 
(m w x m w )-matrix \1/ , (n x x n x )-matrix $, and £u>o structured matrix variables 



S:-- 



' St 
















s 2 _ 








55) 



// the problem l[81\) -(84\) is feasible and the unknown variables have been found, then the output- 
feedback controller gain matrix is determined by K = S^ 1 L 1 . 

The proof is similar to that of [70] where it is derived for the norm performance criterion. 

Proof. Let a solution to the problem (ISTl) - (|gl"l) exist. Performing a congruence transformation 
with 

blockdiag (I mw ,T^,I P J, blockdiag (Tj , I mw , , I Pz ) 
on the inequalities (1821 . (|83|) . respectively, leads to 
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where the plant realization matrices are derived from the backward transformation of (T79|) . Let 
us denote S := T^ST U , $ := Tj$T u . Then from (ISo^) and definition of K = S^Li it follows 
that 
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and the above LMIs can be rewritten as 
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or, in terms of the closed-loop realization (180|) . as 
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Then, performing a congruence transformation with 

blockdiag (I mw , S' 1 , J P2 ) , blockdiag (I n:c , 
on the last inequalities, respectively, we have 
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(86) 



-< 0. 



(87) 



From the inequality 
it is clear that 



Then, by Lemma [U from (18~T1) . ([8~B"j) . f[8"7]) . (I8"~"j) it follows that the controller gain matrix K is 
the solution to Problem [3] for the closed- loop realization (18*01) , which completes the proof. □ 

Remark 2. Unlike the proofs of Theorems[J}^ there is no equivalence between the synthesis 
inequalities l[8l\) -(84\) and the conditions ftlfy-^Lfy of LemmaUl The synthesis LMIs (EI)), 
establich only sufficient conditions for the inequalities < TT3|) . |7j] ) o/ Lemma U\ to be solvable. 
This also concerns a synthesis theorem below. 

Corollary 6. With the notation 7 := r y 2 , the conditions of Theorem^ allow for 7 to be mini- 
mized via solving the convex optimization problem 



minimize 7 
over q?,$,S,L,r],7 satisfying l[81\) -[84\). 



(88) 



If the problem ("~~5J) is solvable, the controller gain matrix is constructed just as in Theorem^ 

Now consider the singular filtering problem when the matrix D yw of the plant (pP) is ze- 
ro and the matrix C y is of full row rank. Then there exists a nonsingular state coordinate 
transformation matrix T y such that [70] 

C y := CyT- 1 = [ I Py 0]. (89) 
Under this transformation, the plant realization matrices become 



A :— T y AT y , 



B„ 



TyB w , 
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TyB u1 



C z '.— C z T y 1 . 



(90) 
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Theorem 5. Suppose that the plant P described by (QP is such that D yw = and rankC^ = p y . 
Given a ^ 0, 7 > 0, a static output-feedback controller defined by ( f50j) solving Problem^ for 
the closed-loop realization 
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exists if the convex problem 
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(92) 
(93) 



^0, 



(94) 



(95) 



where A, B w , C z , C y are defined by (OSJ), is feasible with respect to the scalar variable rj, real 
(m w x m w )-matrix ty, (n x x n x )-matrix IT, and two structured matrix variables 



R :-- 



Rt 
R 2 



M := [ Mi ] 



(96) 



// the problem / fPDj) -/ fP3]) feasible and the unknown variables have been found, then the output- 
feedback controller gain matrix is determined by K = M 1 R^ 1 . 

The proof is dual to that of Theorem H] and similar to that of [70j where it is derived for 
the "Hoc norm performance criterion. 

Proof. Let a solution to the problem (192]) - (195]) exist. Substitute the realization matrices defined 
by (190]) to the LMIs (193]) . (1941) . Perform a congruence transformation with 

blockdiag (I mw , T~, I Pz ) , blockdiag (T^ 1 , J mro , T~ x , J Pa ) 

on the LMIs ([93]), (JS]), respectively. Then define R := T~ 1 RT~ T and IT := T y l TlT y T . From (J96]) 
and definition of K = M^ 1 it follows that MT~ T = KC y R, and the LMIs (l9"3"j), (l9"i]l can be 
rewritten as 
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or, in terms of the closed-loop realization (19T]) . as 
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(97) 



Then, performing a congruence transformation with blockdiag (R T , I mw , I nx , I Pz ) on the last 
inequality leads to 
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(98) 
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From the inequality 
it is clear that 



Let us define $ := IT 1 . Then, by Lemma OQ from $2$, ([97]), Q2B), (JHSD it follows that the 
controller gain matrix K is the solution to Problem [3] for the closed-loop realization ()9ip . which 
completes the proof. □ 

Corollary 7. Wii/i i/ie notation 7 := 7 2 ; £/ie conditions of Theorem^ allow for 7 to be mini- 
mized via solving the convex optimization problem 



minimize 7 
over \1>, n, R, M, 77, 7 satisfying 



(99) 



If the problem IPPj) solvable, the controller gain matrix is constructed just as in Theorem^ 

It is noted in [70] that since the singular control and filtering problems are dual, the convex 
feasibility problems (l8Tj) -( l84|) and (192]) - (195]) of Theorems H] and [5] are in a sense dual too, just as 
the convex optimization problems (188 p and (1991) of Corollaries E] and Replacing the realization 
matrices and the variables under the coordinate transformation in the formulas of Theorem H] 
and Corollary |6] as 



{-^5 Rwi Bill C Zl DzWl Cy, Dy w ^ 

{rj,*,<S>,S,L} 



J/iT fiT (iT dT r)T dT nT 1 



{ij,$,n,i? T ,i T } 



we obtain the respective formulas of Theorem and Corollary [7] and the controller is given by 
K —+K T . 

It is shown in [70] that the results of Theorem H] and Corollary [6] can be applied to synthesis of 
decentralized anisotropic suboptimal and 7-optimal static output-feedback and fixed-order con- 
trollers. In turn, Theorem [5] and Corollary [7] allow to get a solution to simultaneous anisotropic 
output-feedback control problems. These topics are beyond the limits of this paper and may be 
discussed elsewhere. 
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3.6 Fixed-order controller via convex optimization 



It is well-known (sec e.g. [23J) that the fixed-order dynamic controller synthesis problem can 
be embedded into a static output-feedback design problem by augmentation of the plant states 
with the controller states as 
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(100) 



The closed-loop realization is then given by 

1^ [ Cy *Byw ] 

where the gain matrix K incorporates the controller parameters 



A 2 
e £> 



A 



B v 



+ 



B u 



A + B U KCy 

C z -\- D zw KCy 



B w + B U KV 
V zw + V ZU KV 



yw J 



K :-- 



A c B c 

a. D r 



(101) 



Therefore, if the realization of the plant ([T]) has one of the matrices D zu or D yw identically 
zero with B u or C y of full column/row rank, respectively, we can make use of Theorem H] 
and Corollary [6] or Theorem [5] and Corollary [7] to find the fixed-order anisotropic 7-optimal 
(suboptimal) controller as the static output-feedback gain fllOip for the realization (11001) of the 
augmented plant. 



4 Numerical examples 

In this section we provide several purely illustrative numerical examples of the anisotropic 
7-optimal controller design via convex optimization. Only two special design cases are consid- 
ered, namely, the full-order output-feedback controller and static output-feedback gain defined 
in Theorems [2] and El respectively. As regards general Problems [Q, [3] of the anisotropic subop- 
timal controller design with the solutions defined by Corollaries (2J HI testing and benchmark of 
various algorithms for finding reciprocal matrices under convex constraints (e.g., [58, 63, 64J) 
is the issue of future work and will be presented elsewhere. However, it should be mentioned 
that the algorithms of [58| [63] have been tested on some reasonable number of state-space 
realizations randomly generated by the MATLAB Control Systems Toolbox function drss and 
some models from the COMPl e ib collection [7TJ [72]. The numerical experiments have shown 
that application of both of that algorithms often leads to convergence to local minima and 
depend on initial conditions. The randomized technique proposed in [61] aimed at generation 
of the initial conditions seems to be able to improve the situation. 

All computations have been carried out by means of MATLAB 7.9.0 (R2009b), Control 
System Toolbox, and Robust Control Toolbox in combination with the YALMIP interface [55] 
and the SeDuMi solver [54J with CPU P8700 2 x 2.53GHz. 
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4.1 Full-order output-feedback design 



4.1.1 TU-154 aircraft landing 

First we consider the problem of longitudinal flight control in landing approach under the 
influence of both deterministic and stochastic external disturbances in conditions of a windshear 
and noisy measurements. The control aims at disturbance attenuation and stabilization of the 
aircraft longitudinal motion along some desired glidepath. The linearized discrete time-invariant 
model of TU-154 aircraft landing is given in [73], where the problem was solved by means 
of the anisotropic optimal controller derived in [45] . Here we present the results of solving 
the anisotropic 7-optimal full-order synthesis problem via convex optimization as defined in 
Theorem [2] and Corollary |2J 

The mathematical model of the aircraft longitudinal motion defining deviation from a 
nominal trajectory was derived in [73] at the trajectory point characterized by the airspeed 
Vq = 71.375 m/sec, flying path slope angle #0 = —2.7 deg, pitch angle rate u z o = deg/sec, 
pitch angle ^0 = deg, height h = 600 m, and thrust T = 52540 newton. The model has 
order n x = 6, two control inputs (the signal A$ cy generated by the controller to deflect the gen- 
eralized ailerons and the throttle lever position A5 t ) and two measured outputs (the airspeed 
AV + wy,k and the height Ah + Wh,k)- The sampling time of the model At = 0.01 sec. 

The anisotropic 7-optimal controller K a was derived from a solution to the convex opti- 
mization problem (1491) as defined in Theorem [2j The state-space realization of the anisotropic 
7-optimal controller K a computed for the mean anisotropy level a = 0.7 is presented below 
together with the realizations of %2 and optimal controllers K 2 and computed by 
MATLAB Robust Control Toolbox functions h2syn (Riccati equations technique) and hinf syn 
(LMI optimization technique): 



K 2 = 
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The results of simulation of the closed-loop systems in conditions of a windshear and noisy 
measurements are presented together with the problem solution results in Table [T] below and 



23 



Table 1: TU-154 aircraft landing. Comparison of closed-loop systems 







Controller in feedback loop 






K 2 


K a 


Koo 


Solution 


results: 








min7 




0.516 


5.4203 


10.894 


WTzwh 




0.516 


1.1473 


3.1448 


III Tz wl 0.7 




7.8391 


5.1768 


5.5944 


\\Tzw Woo 




15.855 


10.93 


10.891 


CPU time, 


sec 


0.78001 


5.928 


1.7004 


Simulation results: 


max|AV|, 


m/sec 


11.3 


3.559 


4.329 


max | Ah\, 


m 


54.79 


46.87 


39.79 


max \A8 \, 


deg 


14.86 


16.04 


31.6 


max |Ao; z |, 


deg/ sec 


4.884 


5.043 


10.56 


max |Al?|, 


deg 


19.06 


19 


38.08 


max |AT|, 


kN 


7.263 


22.58 


42.48 


max Ai? cy 


deg 


20.7 


20.8 


21.91 


max | ASt | , 


deg 


8.224 


29.25 


29.23 



illustrated in Figures In the simulation we use a typical wind profile described by the ring 
vortex downburst model [74] . 

From the solution results in Table [1] we can conclude that 

• the respective minimum square root values of the objective functions satisfy 72 < 7a < 7c«; 

• the a-anisotropic norm of the closed-loop system with the anisotropic 7-optimal controller 
satisfies |||Tzw|| 0>7 < la] the controller is actually suboptimal. 

Analysis of the simulation results presented in Table [T] and Figures [2HH shows that 

• the anisotropic 7-optimal controller results in the least maximal absolute deviation of the 
airspeed and admissible maximal absolute deviation of the height; 

• the worst maximal absolute deviations of the controlled variables are demonstrated by 
the %2 optimal controller; 

• the anisotropic controller provides the maximal absolute deviation of the thrust required 
for the manoeuvre almost two times less than the additional thrust required by the system 
with the Hoc controller; 

• the same concerns the maximal absolute deviations of the trajectory slope angle, pitch 
rate, and pitch; 

• the least maximal additional thrust is required by the closed-loop system with the %2 
optimal controller; 

• the maximal values of the control signals of the anisotropic and Hoo controllers are close, 
the control generated by the anisotropic controller looks more realistic. 

The anisotropic 7-optimal controller is obviously more effective than the %2 controller and 
less conservative than the "Hco controller in this example of the disturbance attenuation problem. 
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Time (sec) Time {sec] 



Figure 2: TU-154 aircraft landing. Airspeed AV, height Ah (left plots) and control signals 
Ai? C2/ , A5 t (right plots) 




Time (sec) Time (sec) 



Figure 3: TU-154 aircraft landing. Trajectory slope angle A9, pitch angle rate Au z (left plots), 
pitch angle A$, thrust AT (right plots) 
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Figure 4: TU-154 aircraft landing. Wind profile (left plots) and noisy measurements (right 
plots) 

4.1.2 COMPl e ib examples 

The anisotropic 7-optimal full-order controllers have been computed for some models from 
the COMPl e ib collection [7TJ [72] listed below in Table [2J All of them were converted from 
continuous- to discrete-time models with the sampling time At. It is known from [72] that almost 
all of these models (excepting ROC5) are SOF-stabilizable, but here the respective problems 
are solved by the dynamic full-order output-feedback controllers for the testing purpose solely. 
In [53] it is shown that satisfying the conditions of SANBRL with a — > 0, +00 ensures the H2 
and Hoc norms not to exceed a given threshold value. Therefore the I-L2 and controllers for 
the respective problems have also been derived as the limiting cases of the anisotropic controller 
from a solution to the convex optimization problem (149]) as defined in Theorem [2] but with the 
respective input mean anisotropy levels a = and a = +00. 

Below we present the solution and simulation results for the autopilot control problem for 
an air-to-air missile (AC4) initially presented in [6U] , where this problem is considered in robust 
setting and requires that the autopilot generates the tail deflection 5 to produce an angle of 
attack a corresponding to a manoeuvre defined by the guidance law [60]. More precisely, the 
control aims at tracking step input commands a c with a steady state accuracy of 1%, achieving 
a rise time less than 0.2 sec and limiting overshoot to 2% over a range of angles of attack 
±20deg and variations in Mach number 2.5 to 3.5 [60] . The model AC4 from the COMPl e ib 
collection [71] does not take into account the variations in Mach number and therefore not 
include uncertain parameters. 

The state-space realization of the anisotropic 7-optimal controller K a computed for the 
mean anisotropy level a = 0.015 is presented below together with the realizations of the %2 
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Table 2: Examples from the COMPl e ib collection [TTj 172]. Full-order design 
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and "Hoc controllers K 2 and K r 
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The results of simulation of the closed- loop systems with the anisotropic and T^oo 
controllers K 2 , K a and in conditions of noisy measurements are illustrated in Figures 
[HJ In simulation we generated the reference commands a c as steps with random amplitudes 
and equal fixed durations. The step responses in Figure [8] show that the closed-loop rise time 
pointed out in these plots does not exceed the desired 0.2 sec for all three controllers. At that, 
the diagrams of Figure [5] demonstrate acceptable tracking performance and lesser amplitude of 
the control 5 required in the closed-loop system with the anisotropic controller in comparison 
with the "Hoc, controller. 
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Figure 5: Model AC4 (air-to-air missile) [SIl EU E2]- Angle of attack a (top plot) and tail 
deflection angle 8 (bottom plot) 
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Time [sec] 

Figure 6: Model AC4 (air-to-air missile) [HOI EH [T5] . Controlled output z 



28 




a- 1 2 3 4 5 6 7 



Figure 7: Model AC4 (air-to-air missile) [60j EU E2] • Measurement y (top plots) and pitch rate 
sensor noise n (bottom plot) 
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Figure 8: Model AC4 (air-to-air missile) [6Ql ED E2] ■ Bode diagram (top plots) and step response 
(bottom plots) 
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4.2 Static output-feedback design 



The anisotropic 7-optimal static output-feedback controllers have been computed for a number 
of singular filtering problems from the COMPl e ib collection [71] [72] listed below in Table |3j As 
above, all of them were converted from continuous- to discrete-time models with the sampling 
time At. The W2 and 7-Loo controllers for the respective problems have been derived as the lim- 
iting cases of the anisotropic controller from a solution to the convex optimization problem f[9"9"j) 
as defined in Theorem and Corollary [7J but with the respective input mean anisotropy levels 
a = and a = +00. 

Table 3: Examples from the COMPl e ib collection [7TJE2]. Static output-feedback design 



Model 


(n x ,m u ,p y ) 


At 

(sec) 


min 72 
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min 7 a 
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0.01 
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1 
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0.32302 
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(7,2,3) 


0.001 


0.035481 


0.01 
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UWV 


(8,2,2) 


0.001 
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0.03 
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0.024207 
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For the purely illustrative purpose, below we present the solution and simulation results 
for the aircraft control problem (ACl) initially considered in [75]. The model AC1 from the 
COMPlJb collection [TTJ is recast into a disturbance attenuation singular filtering problem 
with noiseless measurements. The anisotropic 7-optimal static gain K a computed for the mean 
anisotropy level a = 0.9 is presented below together with the H2 and %oo gains K 2 and : 



7.278 ■ 10~ J -0.9994 -0.000203 
-0.0002887 -0.002706 -0.9966 
-0.9871 -14.34 49.35 



K a 
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54.56 
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The results of simulation of the closed-loop systems in conditions of a windshear are present- 
ed together with the problem solution results in Table H] below and illustrated in Figures I9TTTT1 
In the simulation we use the same wind profile as in the example of TU-154 aircraft flight 
control in Section 14.1. II 

The solution results presented in Table @] shows that 

• the respective minimum square root values of the objective functions satisfy 72 < 7a < 7oo; 

• the a-anisotropic norm of the closed-loop system with the anisotropic 7-optimal static 
gain satisfies ||Tzw| 0>9 < la', the controller is actually suboptimal; 
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Table 4: Example AC1 (aircraft) from the COMPl e ib collection [7T| 172]. Comparison of closed- 
loop systems 





Controller in feedback loop 




K 2 \ K a \ h\ 



Solution results: 
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Simulation results: 
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max \xi\, 


m 


3.152 


3.412 


3.35 


max \x2\, 


m/sec 


0.1647 


0.1108 


0.124 


max |a;3 1, 


deg 


0.02948 


0.0192 


0.02172 


max |x4|, 


deg/sec 


0.008596 


0.006704 


0.006841 


max |a?s |) 


m/sec 


0.406 


0.278 


0.3097 


max 


10" 1 deg 


0.1648 


0.1108 


0.124 


max |«2|i 


m/sec 2 


0.0299 


0.01922 


0.02173 


max |«3|, 


deg 


0.2117 


0.1355 


0.154 



• the "H 2 and Hoc norms of the closed-loop systems with the respective 7-optimal gains 
satisfy ||7zvk||2 < 72, ||7zjy||oo < 7oo! the %2 and controllers are actually suboptimal 
too. 

The simulation results presented in Table @] and Figures IUHTT1 allow to conclude that 

• the anisotropic 7-optimal output-feedback static gain leads to the least maximal absolute 
deviations of the forward speed x 2 , pitch angle x 3 , pitch angle rate x 4 , and vertical speed 
X5, at that the least maximal absolute deviation of the height error x\ is achieved with 
the %2 7-optimal static gain; 

• the worst maximal absolute values of the controlled output are demonstrated by %2 7~ 
optimal static gain; 

• the anisotropic 7-optimal static gain leads to the least maximum absolute amplitudes of 
the control signals. 

5 Conclusion 

In this paper, we have proposed a solution to the anisotropic suboptimal and 7-optimal con- 
troller synthesis problems by convex optimization technique. The anisotropic suboptimal con- 
troller design is a natural extension of the optimal approach developed in [15]. Instead of 
minimizing the anisotropic norm of the closed-loop system, the suboptimal controller is only 
required to keep it below a given threshold value. The general fixed-order synthesis procedure 
employs solving an inequality on the determinant of a positive definite matrix and two linear 
matrix inequalities in reciprocal matrices which make the general optimization problem non- 
convex. By applying the known standard convexification procedures it have been shown that 
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Figure 9: Model AC1 (aircraft) [75j [7U [72]. Controlled output and wind profile (left plots), 
control (right plots) 
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Figure 10: Model AC1 (aircraft) [75], EH E2] . Height error x\, forward speed 22, vertical speed 
25 (left plots), pitch angle 23, pitch angle rate 24 (right plots) 
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the resulting optimization problem can be made convex for the full-information state-feedback, 
output-feedback full-order controllers, and static output-feedback controller for some specific 
classes of plants defined by certain structural properties. In the convex cases, the anisotropic 
7-optimal controllers are obtained by minimizing the squared norm threshold value subject 
to convex constraints. In comparison with the solution to the anisotropic optimal controller 
synthesis problem derived in [45] which results in a unique full-order estimator-based controller 
defined by a complex system of cross-coupled nonlinear matrix algebraic equations, the proposed 
optimization-based approach is novel and does not require developing specific homotopy-like 
computational algorithms. 



References 

[1] J.C.Doyle. Guaranteed margins for LQG regulators. IEEE Trans. AC, 23:756-757, 1978. 

[2] J.C.Doyle, K. Glover, P.P. Khar gonekar, and B.A.Francis. State-space solutions to stan- 
dard H 2 and control problems. IEEE Trans. AC, 34:831-848, 1989. 

[3] D.H. Jacobson. Extensions of Linear- Quadratic Control, Optimization and Matrix Theory. 
Academic Press, NY, 1977. 

[4] P. Whittle. Risk-sensitive linear/quadratic/Gaussian control. Adv. Appl. Prob., 13:764- 
777, 1981. 

[5] P. Whittle. Entropy-minimizing and risk-sensitive control rules. Syst. & Contr. Lett., 13: 
1-7, 1989. 



33 



[6] D.S. Bernstein and W.M. Haddad. LQG control with an Hoc performance bound: a Riccati 
equation approach. IEEE Trans. AC, 34:293-305, 1989. 

[7] D. Mustafa and K. Glover. Minimum Entropy Hoc Control. Springer- Verlag, NY, 1991. 

[8] K. Glover and J.C. Doyle. State-space formulae for all stabilizing controllers that satisfy an 
"Hoo-norm bound and relations to risk sensitivity. Syst. & Contr. Lett., 11:167-172, 1988. 

[9] P.A. Iglesias, D. Mustafa, and K. Glover. Discrete time V,^ controllers satisfying a mini- 
mum entropy criterion. Syst. & Contr. Lett, 14:275-286, 1990. 

[10] D. Mustafa, K. Glover, and D. Limebeer. Solutions to the Hoo general distance problem 
which minimize an entropy integral. Automatica, 27:193-199, 1991. 

[11] P.A. Iglesias and D.Mustafa. State-space solution of the discrete-time minimum entropy 
control problem via separation. IEEE Trans. AC, 38:1525-1530, 1993. 

[12] I. Yaesh and U. Shaked. Minimum entropy static output-feedback control with an 
^oo-norm performance bound. IEEE Trans. AC, 42:853-858, 1997. 

[13] E. Fridman and U. Shaked. Robust minimum entropy static output-feedback control 
of singularly perturbed systems. Automatica, 36:1181-1188, 2000. 

[14] K. Zhou, K. Glover, B.A. Bodenheimer, and J.C. Doyle. Mixed %2 and "H^ performance 
objectives I: Robust performance analysis, II: Optimal control. IEEE Trans. AC, 39: 
1564-1574, 1575-1587, 1994. 

[15] H. Rotstein and M. Sznaier. An exact solution to general four-block discrete-time mixed 
^/Uoo problems via convex optimization. IEEE Trans. AC, 43:1475-1481, 1998. 

[16] R. Miradore and G. Ricci. Mixed ^2/^00 control: the discrete-time case. Syst. & Contr. 
Lett, 54:1-13, 2005. 

[17] P.P. Khargonekar and M.A. Rotea. Mixed T^/'Hoo control: a convex optimization approach. 
IEEE Trans. AC, 36:824-837, 1991. 

[18] S.P. Boyd and L.E1 Ghaoui. Method of centers for minimizing generalized eigenvalues. Lin. 
Alg. AppL, 188:63-111, 1993. 

[19] Yu. Nesterov and A. Nemirovskii. Interior-point Polynomial Methods in Convex Program- 
ming, volume 13 of Studies in Applied Mathematics. SIAM, Philadelphia, PA, 1994. 

[20] A. Nemirovskii and P. Gahinet. The projective method for solving linear matrix inequali- 
ties. Math. Programming Series B, 77:163-190, 1997. 

[21] S. Boyd, L. El Ghaoui, E. Feron, and V. Balakrishnan. Linear Matrix Inequalities in Sys- 
tems and Control theory. SIAM, Philadelphia, PA, 1994. 

[22] P. Gahinet and P. Apkarian. A linear matrix inequality approach to "Hoc control. Int. J. 
of Robust and Nonlinear Contr., 4:421-448, 1994. 



34 



[23] T. Iwasaki, R.E. Skelton. All controllers for the general "Hoc control problem: LMI existence 
conditions and state space formulas. Automatica, 30:1307-1317, 1994. 

[24] C.W. Scherer. Multiobjective % 2 /^°o control. IEEE Trans. AC, 40:1054-1062, 1995. 

[25] X. Chen and J.T. Wen. A linear matrix inequality approach to the general mixed T^/^oo 
control problem. Proc. American Control Conf., pages 1443-1447, 1995. 

[26] C.W. Scherer, P. Gahinet, and M. Chilali. Multiobjective output-feedback control via LMI 
optimization. IEEE Trans. AC, 42:896-911, 1997. 

[27] I. Masubuchi, A. Ohara, and N. Suda. LMI-based controller synthesis: A unified formulation 
and solution. Int. J. of Robust and Nonlinear Contr., 8:669-686, 1998. 

[28] H.A.Hindi, B.Hassibi, and S.P.Boyd. Multiobjective %2/%00-optimal control via finite 
dimensional Q-parametrization and linear matrix inequalities. Proc. American Control 
Conf., pages 3244-3248, 1998. 

[29] M.C. Oliveira, J.C. Geromel, and J. Berbnussou. An LMI optimization approach to multi- 
objective controller design for discrete-time systems. Proc. IEEE Conf. on Decision and 
Control, pages 27-38, 1999. 

[30] C.W. Scherer. An efficient solution to multi-objective control problems with LMI objec- 
tives. Syst. & Contr. Lett, 40:43-57, 2000. 

[31] P. Apkarian, P.C. Pellanda, and H.D. Tuan. Mixed T^/^oo multi-channel linear parameter- 
varying control in discrete time. Syst. & Contr. Lett., 41:333-346, 2000. 

[32] C.W. Scherer. Multi-objective control without Youla parametrization. In Perspectives in 
Robust Control. Lecture Notes on Control and Information Sciences, 268:311-325, 2001. 

[33] D.Arzelier and D.Peaucelle. An iterative method for mixed % 2 /%oo synthesis via static 
output feedback. Proc. IEEE Conf. Dec. Contr., pages 3464-3469, 2002. 

[34] J. Yu. A new static output feedback approach to the suboptimal mixed ^2/^00 problem. 
Int. J. of Robust and Nonlinear Contr., 14:1023-1034, 2004. 

[35] S.V. Gusev. Minimax control under a bound on the partial covariance sequence of the 
disturbance. Automatica, 31:1287-1301, 1995. 

[36] S.V. Gusev. Minimax control under a restriction on the moments of disturbance. Proc. 
34th IEEE Conf. on Decision and Control, New Orleans, USA, pages 1195-1200, 1995. 

[37] S.V. Gusev. Method of moment restrictions in robust control and filtering. Proc. 13th 
IFAC World Congress, San-Francisco, USA, pages 415-420, 1996. 

[38] C.W. Scherer. Robust controller design by output feedback against uncertain stochastic 
disturbances. Proc. 3rd IFAC Symp. on Robust Control Design, Prague, Czechia, 2000. 



35 



[39] A.V. Semyonov, I.G. Vladimirov, and A. P. Kurdjukov. Stochastic approach to Hoo- 
optimization. Proc. 33rd IEEE Conf. on Decision and Control, Florida, USA, pages 2249- 
2250, 1994. 

[40] I.G. Vladimirov, A. P. Kurdjukov, and A.V. Semyonov. Anisotropy of signals and the en- 
tropy of linear stationary systems. Doklady Math., 51:388-390, 1995. 

[41] I.G. Vladimirov, A. P. Kurdjukov, and A.V. Semyonov. The stochastic problem of Hoo- 
optimization. Doklady Math., 52:155-157, 1995. 

[42] I.G. Vladimirov, A. P. Kurdjukov, and A.V. Semyonov. On computing the anisotropic norm 
of linear discrete-time-invariant systems. Proc. 13th IFAC World Congress, San-Francisco, 
USA, pages 179-184, 1996. 

[43] P. Diamond, I.G. Vladimirov, A. P. Kurdyukov, and A.V. Semyonov. Anisotropy-based per- 
formance analysis of linear discrete time invariant control systems. Int. J. of Contr., 74: 
28-42, 2001. 

[44] I.G. Vladimirov, P. Diamond, and P. Kloeden. Anisotropy-based performance analysis of 
finite horizon linear discrete time varying systems. Automat. & Remote Contr., 8:1265- 
1282, 2006. 

[45] I.G. Vladimirov, A. P. Kurdjukov, and A.V. Semyonov. State-space solution to anisotropy- 
based stochastic "Hoc-optimization problem. Proc. 13th IFAC World Congress, San Fran- 
cisco, USA, pages 427-432, 1996. 

[46] A. P. Kurdyukov and E.A. Maximov. State-space solution to stochastic Hoo-optimization 
problem with uncertainty. Proc. 16th IFAC World Congr., Prague, Czechia, 2005. 

[47] A. P. Kurdyukov and E.A. Maximov. Solution of the stochastic "Hoo-optimization problem 
for discrete time linear systems under parametric uncertainty. Automat. & Remote Contr., 
67:1283-1310, 2006. 

[48] P. Diamond, A. P. Kurdjukov, A.V. Semyonov, and I.G. Vladimirov. Homotopy methods 
and anisotropy-based stochastic Hoo optimization of control systems. Report 97-14 of The 
University of Queensland, Australia, pages 1-22, 1997. 

[49] CD. Charalambous and F. Rezaei. Stochastic uncertain systems subject to relative entropy 
constraints: induced norms and monotonicity properties of minimax games. IEEE Trans. 
AC, 52:647-663, 2007. 

[50] A. P. Kurdyukov, E.A. Maximov, and M.M. Tchaikovsky. Anisotropy-based bounded real 
lemma. Proc. 19th Int. Symp. on Mathematical Theory of Networks and Systems, Bu- 
dapest, Hungary, 2010, pages 2391-2397. 

[51] E.A. Maximov, A. P. Kurdyukov, and I.G. Vladimirov. Anisotropic norm bounded real lem- 
ma for linear discrete time varying systems. Proc. 18th IFAC World Congr., Milano, Italy, 
2011. 



36 



[52] M.M. Tchaikovsky, A. P. Kurdyukov, and V.N. Timin. Strict anisotropic norm bounded real 
lemma in terms of inequalities. Proc. 18th IFAC World Congr., Milano, Italy, 2011. 

[53] M.M. Tchaikovsky, A. P. Kurdyukov, and V.N. Timin. A convex formulation 

of strict anisotropic norm bounded real lemma. Preprint. Available from 
http : //arxiv . org/abs/1108 . 5140, 2011. 

[54] J.F.Sturm. Using SeDuMi 1.02, a MATLAB toolbox for optimization over symmetric 
cones. Optimization Methods and Software, 11:625-653, 1999. 

[55] J. Lofberg. YALMIP: A toolbox for modeling and optimization in MAT- 
LAB. Proc. of the CACSD Conference, Taipei, Taiwan, 2004. Available from 
control . ee . ethz . ch/~joloef /wiki/pmwiki .php. 

[56] A. Ben-Tal and A. Nemirovskii. Lectures on Modern Convex Optimization. Technion, Haifa, 
Israel, 2000. 

[57] T. Iwasaki and R.E. Skelton. The W-centering algorithm for the dual LMI Problem: A 
new approach to fixed order design. Int. J. of Contr., 62:1257-1272, 1995. 

[58] L. El Ghaoui, F. Oustry, and M.A. Rami. A cone complementary linearization algorithm 
for static output-feedback and related problems. IEEE Trans. AC, 42:1171-1176, 1997. 

[59] P. Apkarian and H.D. Tuan. Concave programming in control theory. J. of Glob. Opt., 15: 
343-370, 1999. 

[60] B. Fares, P. Apkarian, and D. Noll. An augmented Lagrangian method for a class of LMI- 
constrained problems in robust control theory. Int. J. of Contr., 74:348-360, 2001. 

[61] P. Apkarian, D. Noll, and H.D. Tuan. Fixed-order 'H OQ control design via a partially aug- 
mented Lagrangian method. Int. J. of Nonlinear and Robust Contr., 13:1137-1148, 2003. 

[62] D.Noll, M. Torki, and P. Apkarian. Partially augmented Lagrangian method for matrix 
inequality constraints. SI AM J. on Opt, 15:161-184, 2004. 

[63] D.V. Balandin and M.M. Kogan. Synthesis of controllers on the basis of a solution of linear 
matrix inequalities and a search algorithm for reciprocal matrices. Automat. & Remote 
Contr., 66:74-91, 2005. 

[64] B.T. Polyak and E.N. Gryazina. Hit-and-Run: Randomized technique for control problems 
recasted as concave programming. Proc. 18th IFAC World Congr., Milano, Italy, 2011. 

[65] B.T. Polyak and E.N. Gryazina. Markov chain Monte Carlo method exploiting barrier func- 
tions with applications to control and optimization. Proc. IEEE Multi-Conf. on Systems 
and Control, pages 1553-1557, 2010. 

[66] P. Gahinet. Explicit controller formulas for LMI-based T-ioo synthesis. Automatica, 32: 
1007-1014, 1996. 

[67] A.S. Poznyak. Advanced Mathematical Tools for Automatic Control Engineers. Volumes 
1,2: Deterministic Techniques, Stochastic Techniques. Elsevier, 2008, 2009. 



37 



[68] V. Kucera. Stability of discrete linear feedback systems. Proc. IFAC World Congress, 
Boston, Massachusetts, paper No. 44-1, 1975. 

[69] D.C.Youla, H.A. Jabr, and J. J. Bongiorno. Modern Wiener-Hopf design of optimal con- 
trollers. Part II: the Multivariable case. IEEE Trans. AC, 21:319-338, 1976. 

[70] K.H. Lee, J.H. Lee, and W.H. Kwon. Sufficient LMI conditions for T-Loo output feedback 
stabilization of linear discrete-time systems. IEEE Trans. AC, 51:675-680, 2006. 

[71] F. Leibfritz. COMPl e ib: Constraint matrix-optimization problem library — a collection 
of test examples for nonlinear semidefinite programs, control system design and related 
problems. Tech. rep. of the University of Trier, Germany, 2004, http://www.complib.de. 

[72] F. Leibfritz and W. Lipinski. Description of the benchmark examples in COMPl e ib 1.0. 
Tech. rep. of the University of Trier, Germany, 2003, http://www.complib.de. 

[73] A. P. Kurdyukov, B.V.Pavlov, V.N. Timin, and I.G. Vladimirov. Longitudinal anisotropy- 
based flight control in a wind shear. Proc. 16th IFAC Symposium on Automatic Control 
in Aerospace, Saint-Petersburg, Russia, 2004. 

[74] M. Ivan. A ring vortex downburst model for flight simulation. J. Aircraft, 23:232-236, 
1996. 

[75] Y.S. Hung and A.G.J. MacFarlane. Multivariable feedback: A quasi- classical approach, vol- 
ume 40 of Lecture Notes in Control and Information Sciences. Springer- Verlag, Berlin, 
Heidelberg, New York, 1982. 



38 



